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ABSTRACT 


We  clearly  elucidate  the  kinematic  constraints,  and  the  corresponding  kinematic 
indeterminacy  of  part  of  the  deviatoric  stress  tensor,  in  a  rigid-viscoplastic  single  crys¬ 
tal  lacking  five  independent  slip  systems.  The  indeterminate  stress  component  is  a 
Lagrange  multiplier  enforcing  the  kinematic  constraint,  and  it  must  be  determined 
from  equilibrium  considerations.  A  simple  polycrystalline  model  is  constructed  which 
precisely  satisfies  iocal  kinematic  constraints  as  well  as  global  compatibility.  Volume- 
average  global  stresses  are  obtained  by  approximating  the  local  constraint  stress  as 
the  corresponding  projection  of  the  (to-be-determined)  global  stress.  Applications  of 
the  model  to  hexagonal  crystals  without  pyramidal  slip,  and  to  large  deformation  and 
texturing  of  orthrhombic  polycrystalline  materials  (olivine;  HDPE)  are  made. 

l.  Introduction 

PHYSICALLY-BASED  continuum  modeling  of  plastic  deformation  in  ^olycrystalline  ag¬ 
gregates  is  a  broad  and  active  field  of  research.  The  simplest  eff<  e  models  of  this 
sort  are  based  on  the  suggestion  of  TAYLOR  (1938)  that  local  detormation  within  a 
heterogeneous  polycrystalline  aggregate  is  often  approximately  homogeneous  when  a 
representative  material  element  is  subject  to  loading  consistent  with  macroscopically 
uniform  deformation.  This  assumption  assures  both  local  and  global  compatibility  but 
in  general  violates  local  traction  equilibrium  across  grain  boundaries.  In  many  cases, 
however,  major  features  of  both  the  macroscopic  constitutive  behavior  of  the  aggregate 
and  the  evr  ttion  of  crystallographic  texture  with  on-going  deformation  can  be  accu¬ 
rately  s'Tr.i  M  by  TAYLOR-type  models.  When  applied  to  crystalline  deformation 
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by  slip,  various  adaptations  of  the  TAYLOR  model  have  provided  assessments  of  effects 
of  differing  deformation  resistances  for  different  classes  of  slip  systems  (HUTCHINSON, 
1976;  1977),  of  self  and  latent  hardening  of  slip  systems  (ASARO  and  NEEDLEMAN, 
1985),  of  strain  rate  sensitivity  of  slip  (HUTCHINSON,  1976;  1977),  and  of  finite  defor¬ 
mation  elastic/viscoplastic  flow  (ASARO  and  NEEDLEMAN,  1985).  Among  the  effects 
which  cannot  be  directly  addressed  by  such  models  are  the  influences  of  morphological 
texture,  including  grain  shape  and  orientation. 

The  greatest  successes  of  TAYLOR  models  have  occurred  in  cubic  metal  crystals, 
where  the  large  number  of  equivalent  slip  systems  which  (apart  from  modest  differences 
due  to  differential  strain  hardening)  possess  equivalent  deformation  resistances.  The 
flow  envelope  of  a  crystal  can  be  defined  as  the  locus  of  (deviatoric)  stress  components 
required  to  generate  resolved  shear  stresses  equal  to  the  deformation  resistances  on  a 
sufficient  set  of  slip  systems  so  as  to  accommodate  arbitrary  deformation  rates.  For  high 
symmetry  cubic  metal  crystals,  the  flow  envelope  is  roughly  isotropic.  In  crystals  of 
lower  symmetry,  where  the  deformation  resistances  of  different  classes  of  slip  systems 
can  vary  markedly,  the  stress  required  to  activate  slip  systems  accommodating  an 
arbitrary  deformation  rate  can  become  exceedingly  high  for  certain  orientations,  and 
the  flow  envelope  can  become  highly  distorted.  For  this  reason,  the  corresponding 
polycrystalline  aggregate  stress  levels  obtained  from  TAYLOR  type  models  also  become 
large.  In  the  limiting  case  of  unbounded  deformation  resistance  on  all  but  four  or  fewer 
linearly  independent  slip  systems,  the  TAYLOR  estimate  of  aggregate  stress  generally 
becomes  unbounded. 

HUTCHINSON  (1977)  noted  that  the  self-consistent  method  of  HILL  (1965),  as  im¬ 
plemented  by  HUTCHINSON  (1970;  1976)  could  provide  reasonable  estimates  of  macro- 
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scopic  flow  strength  for  aggregates  of  crystals  possessing  slip  systems  of  widely  differing 
deformation  resistance,  including  the  limiting  case  of  less  than  five  linearly  independent 
slip  systems.  The  main  reason  that  the  self-consistent  model  can  deal  with  these  cases 
is  that  it  imposes  much  less  stringent  kinematical  conditions  on  the  deformation  of 
individual  crystals  than  does  the  TAYLOR  model,  with  favorably  oriented  crystals  dis¬ 
playing  higher  deformation  than  average  and  unfavorably  oriented  crystals  displaying 
less  deformation  than  average.  HUTCHINSON  (1977)  applied  the  self-consistent  method 
to  isotropic  distributions  of  hexagonal  crystals.  Recently  TAKESHITA,  et  a/.,  (1989),  us¬ 
ing  the  self-consistent  viscoplastic  approach  of  MOLINARI,  et  at.,  (1987),  modeled  the 
development  of  crystallographic  texture  in  olivine,  an  orthorhombic  mineral  possessing 
only  three  independent  slip  systems. 

While  many  crystalline  materials  display  less  than  five  independent  deformation 
mechanisms  (see,  for  example,  GROVES  and  KELLEY,  1963),  our  interests  in  kinemati¬ 
cally  deficient  crystals  are  motivated  by  semi-crystalline  polymers  such  as  high  density 
polyethylene  (HDPE),  Nylon-6,  and  polyethylene  teraphthalate  (PET).  The  crystalline 
domains  of  these  materials,  while  belonging  to  different  crystal  classes  (orthorhombic, 
monoclinic,  and  triclinic,  respectively),  nonetheless  share  the  feature  that  material  el¬ 
ements  aligned  with  the  lattice  directions  corresponding  to  the  stiff  covalent  bonding 
of  the  backbone  polymer  chain  are  plastically  inextensible.  That  is,  all  known  mecha¬ 
nisms  of  plastic  deformation  in  these  materials,  including  slip  (and,  where  applicable, 
twinning  and  martensitic  transformation),  leave  material  distances  along  the  chain  axis 
invariant.  Large  plastic  deformation  of  semi-crystalline  polymers  can  impart  significant 
crystallographic  texture  to  these  materials.  A  common  textural  feature  is  that  the  crys¬ 
tallographic  chain  axes  preferentially  align  with  the  material  axes  of  maximum  stretch 
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(BOWDEN  and  Young,  1974;  HAUDIN,  1982).  As  this  process  occurs,  the  macroscopic 
stress  required  to  sustain  plastic  flow  can  ultimately  rise  sharply  as  the  locally  available 
sites  for  continued  macroscopic  extension  are  geometrically  deactivated. 

Our  attempts  to  model  large  deformation  in  HDPE  using  the  self-consistent  formu¬ 
lation  of  MOLINARI,  et  al.,  (1987)  provided  qualitatively  correct  deformation  textures, 
but  unsatisfactory  macroscopic  stress-strain  behavior;  macroscopic  average  deviatoric 
stress  components  showed  extreme  fluctuations,  especially  at  sharper  textures.  Be¬ 
cause  of  the  local  inextensibility  constraint,  the  incompressible  viscoplastic  formula¬ 
tion  of  MOLINARI,  et  al.,  (1987)  cannot  fully  resolve  the  local  deviatoric  stress  within 
a  crystallographic  inclusion,  and  the  erratic  macroscopic  stresses  reflect  this  local  inde¬ 
terminacy.  It  seems  that  HUTCHINSON  (1977)  avoided  this  difficulty  by  using  a  slightly 
different  numerical  formulation  of  the  self-consistent  model,  including  a  slight  elastic 
compressibility. 

Here  we  clearly  identify  the  kinematically  indeterminate  component  of  stress  de- 
viator  in  a  rigid-viscoplastic  crystal  possessing  a  kinematic  rank  deficiency  ( e.g .,  a 
locally  inextensible  direction).  This  portion  of  local  stress  state  must  be  determined 
from  equilibrium  considerations,  and  procedures  for  so  doing  are  similar  to  (but  not 
identical  with)  those  for  obtaining  hydrostatic  pressure  in  incompressible  materials. The 
focus  of  the  paper  is  the  development  of  a  new  hybrid  formulation  to  estimate  the  local 
deviatoric  constraint  stresses  in  terms  of  the  (to-be-determined)  macroscopic  stress. 
The  complete  system  of  equations,  which  we  term  the  constrained  hybrid  (CH)  model, 
possesses  a  simple  mathematical  structure.  The  model  is  closely  akin  to  the  simple 
TAYLOR  model  for  crystals  of  full  kinematic  rank  so  that,  for  example,  rate  problems 
of  imposed  macroscopic  deformation  can  be  solved  without  global  iterations,  in  con- 
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trast  to  more  elaborate  self-consistent  models.  Large  deformation  texture  predictions 
and  macroscopic  stress-strain  behavior  are  readily  obtained. 

We  compare  the  model  with  calculations  of  HUTCHINSON  (1977)  for  isotropic  tex¬ 
tures  of  hexagonai  crystals  possessing  only  four  independent  systems.  We  generalize  the 
CH  model  to  cases  of  three  or  fewer  slip  systems  (two  or  more  deviatoric  constraints) 
and  compare  results  of  the  model  with  the  texture  predictions  of  TAKESHITA,  et  al., 
(1989)  in  olivine.  Finally,  we  apply  the  CH  model  in  a  highly  idealized  way  to  the 
deformation  of  initially  isotropic  HDPE. 

Notation  is  based  on  the  following  conventions.  Scalars  are  in  mathematical  italics 
(A,  a ,  a),  vectors  are  lower  case  bold-face  (a),  second  order  tensors  are  upper  case 
bold-face  (A),  and  fourth  order  tensors  are  in  upper  case  caligraphic  type  (A).  Tensor 
(dyadic)  products  are  indicated  by  ‘®’  and  tensor  scalar  products  of  appropriate  order 
are  denoted  by  a  raised  dot.  The  superscripts  uf”,uTn,  and  “-1”  denote  “deviatoric 
part  of”,  “transpose”,  and  “inverse”,  respectively.  When  required,  repeated  cartesian 
subscripts  are  summed  from  1  to  3.  For  example,  on  orthonormal  basis  vectors  e<  : 


c  =  c.e, 
a  =  o,e, 


=»  a  •  c  =  Oi  Cj 


A  —  G,  0  6j 

B  =  B,;e,  ®  e, 


Ac  —  AijCj  ©,• 

AB  —  AikBkj  6j  ®  6j ■ 

A  •  B  =  tr AT5t  =  AijBij 


A  =  B  ®  C  f  A  [G]  =  B  (C  •  G) 

B=D®E  ^  \  AB  =  (C-D)B®E. 
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2.  Constrained  Single  Crystal 


Consider  a  crystal  possessing  N  slip  systems,  of  which  the  ath  possesses  slip  plane 
normal  na  and  slip  direction  s°,  with  na  •  sa  =  0.  Let  the  shear  strain  rate  on  system 
a  be  7®.  In  anticipation  of  large  deformation  we  neglect  elasticity  and  take  D,  the 
traceless  rate  of  deformation  generated  by  these  slips,  to  be  given  by 

N  1 

D  =  Hi*  ~{na®8a  +  aa®  na} 

a=l  ^ 

s  Ei-R”,  (i) 

a=l 

where  R“  is  the  symmetric,  traceless  Schmid  tensor  associated  with  slip  system  a. 

We  further  assume  that  the  span  of  the  N  tensors  R“  has,  in  symmetric  deviatoric 
second  order  tensor  space,  T,,  a  kinematic  deficiency  of  rank  one;  that  is,  a  single 
constraint  (in  addition  to  incompressibility)  exists  on  D.  As  a  physical  example  of 
such  a  constraint,  suppose  that  material  fibers  parallel  to  unit  vector  c  are  inextensible. 
Such  a  constraint  can  arise  in  hexagonal  crystals  when  pyramidal  slip  is  absent  and  only 
basal  and  prismatic  slip  are  available;  in  this  case  the  inextensible  fiber  coincides  with 
[0001]  directions,  HUTCHINSON  (1977).  The  inextensibility  constraint  can  be  written 
as 

c  Dc  =  0.  (2) 

In  order  for  this  to  be  so  for  arbitrary  shears  ia,  it  is  necessary  that  c  -R“c  =  0  for  each 
system.  These  constraints  can  also  be  phrased  in  terms  of  the  tensor  scalar  product 
of  D  (or  R“)  with  the  tensorial  internal  variable  C  =  c  ®  c.  The  tensor  C  is  not 
irreducible,  and  we  find  it  more  convenient  to  retain  the  traceless  symmetric  second 
order  tensor  C'  =  C  -  |I,  where  I  is  the  second  order  identity  tensor,  as  local  internal 
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variable.  Because  D  and  R“  are  traceless,  the  inextensibility  constraint  (2)  can  also 
be  expressed  as 


implying 


D  •  C'  =  0, 

R“  C'  =  0. 


(3) 

(4) 


In  view  of  the  constraints  (4),  we  will  denote  the  span  of  Ra  by  TjXC,  a  four  dimensional 
subspace  of  Tj  which  is  orthogonal  to  the  one  dimensional  subspace  spanned  by  C'. 


We  now  specify  viscoplastic  constitutive  equations  for  the  shears,  7“.  Following 
many  others  ( e.g .,  HUTCHINSON,  1976;  PAN  and  RICE,  1983;  ASARO  and  NEEDLE- 
MAN,  1985),  we  assume  that  7“  depends  on  the  resolved  shear  stress  through  a  power 
law  relation  of  the  form 


Ta  ra  n-1 
•  at  •  *  1*1 

1  ’ 


(5) 


where  n  is  the  rate  exponent,  ga  is  the  slip  system  hardness  (or  shear  strength),  70  is 
a  pre-multiplying  strain  rate,  and  r“  (<  <?“),  the  resolved  shear  stress  on  system  a,  is 
given  by 


ra  =  T  Ra, 


(6) 


where  T  is  the  Cauchy  stress.  Since  R“  is  deviatoric,  we  can  also  express  the  resolved 
shear  stress  as  r“  =  S  •  Ra,  where  the  stress  deviator  is  S  =  X  —  jtr(T)  I.  Slip  system 
strength  may  evolve  with  deformation  according  to 

y°  =  X>a*  1^1,  (7) 

where  ha 0  is  a  matrix  of  hardening  moduli. 
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The  stress  deviator  S  can  be  uniquely  decomposed  els 

S  =  S*  +  ScC' 

=  S-  +  ?  (S  •  C')  C',  (8) 

where  S’  €  TjXC;  that  is,  S*  •  C'  =  0.  Using  (6)  and  (8),  the  resolved  shear  stress  can 
also  be  expressed  as 

r*  =  S-R«  +  ?(SC')C'R* 

=  S*  •  R“,  (10) 

because  of  (4).  When  (10)  is  inserted  into  (5),  and  the  expression  for  7“  inserted  into 
(1),  the  deformation  rate  can  be  written  as 

d  = 

=  X(S-)(S-|.  (11) 

The  deformation  rate  is  independent  of  the  part  of  the  stress  deviator  lying  outside 
Tjxc-  Many  models  of  crystallographic  slip  neglect  effects  of  hydrostatic  stress  on  slip 
rates.  Thus  the  independence  of  the  deformation  rate  from  a  part  of  the  stress  tensor 
is  not  entirely  without  precedent.  In  the  formulation  of  boundary  value  problems  of 
incompressible  crystalline  slip,  the  incompressibility  condition  I  •  D  =  0  is  viewed  as  a 
constraint  and  the  hydrostatic  stress  o  =  I  •  T/  (I  •  I)  =  ^  I  •  T  is  a  Lagrange  multiplier 
enforcing  the  constraint,  the  latter  being  determined  from  equilibrium  considerations. 
Analogously,  in  kinematically  deficient  crystals,  D  •  C1  =  0  is  a  constraint,  and  the 
enforcing  Lagrange  multplier  Sc  —  C  -T/  (Cf  •  C')  =  |  C*  -S  must  be  determined  from 
equilibrium. 
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Thus,  there  are  two  main  features  which  distinguish  the  viscoplastic  response  of  the 
constrained  crystal  from  that  of  “traditional”  crystals.  First,  the  deformation  rate  lies 
in  a  four-dimensional  subspace,  of  deviatoric  space,  T *,  rather  than,  in  gene~al, 

spanning  the  latter.  Secondly,  the  deformation  does  not  depend  on  the  entire  deviatoric 
stress,  but  rather  on  only  the  part  S*  lying  in  a  four-dimensional  deviatoric  subspace. 
Given  these  facts,  it  is  clear  that  the  convexity  (implying  D/S  invertibility)  and  shear 
rate  uniqueness  properties  (HILL,  1956)  usually  associated  with  power  law  viscoplas¬ 
tic  single  crystal  slip  models  («.£.,  HUTCHINSON,  1976)  fail  to  hold.  However,  local 
convexity  connecting  the  respective  four-dimensional  subspaces,  implying  D/S*  invert¬ 
ibility,  does  obtain.  We  will  now  use  these  properties  and  observations  regarding  the 
constrained  single  crystal  model  in  formulating  a  polycrystalline  model  of  deformation. 


3.  POLYCRYSTALLINB  AGGREGATE 


S.l  Kinematics 

Consider  an  aggregate  of  crystals  of  the  type  described  in  the  previous  section  which 
occupy  the  domain  fl  of  volume  Vh.  We  imagine  that  for  points  x  on  the  boundary 
dft,  the  velocity  v  is  prescribed  according  to 

v  =  Lx  =  (D  +  W)  x.  (12) 

Here  L  is  a  constant  macroscopic  “velocity  gradient”  tensor  with  symmetric  part 
D  =  ^  +  LT)  and  skew  part  W  =  Ht+tT)-  Ir  compressibility  dictates  that 

trlL  =  trt)  =  0. 


Within  an  individual  crystal,  let  the  axis  of  inextensibility  be  parallel  to  the  unit 
vector  c.  The  local  deformation  rate  D  within  such  a  crystal  must  satisfy  the  constraints 
(3,4). 


Following  the  insight  of  TAYLOR,  we  would  like  to  choose  a  local  D  which,  while 
satisfying  (3,4),  is  nonetheless  “close”  to  the  macroscopically  imposed  D.  A  minimal 
requirement  of  global  compatibility  within  the  aggregate  ( c.g HILL,  1972)  is  that  the 
volume  average  deformation  rate,  given  by 

<  D  >  =  —  f  Dd7,  (13) 

Vn  Jv n 


equals  D: 


<  D  >  =  D. 


(14) 
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In  order  to  choose  local  values  of  D  which  satisfy  (3,14),  let  us  first  introduce  the 
“intermediate”  traceless  deformation  rate  D,  given  by 

D  =  K  [D] ,  (15) 

where  K  is  a  constant  fourth  order  tensor  (to  be  specified  below)  mapping  Tj  into  itself. 


We  propose  to  express  D  as 

D  =  D  -  <f  C', 

where  d  is  a  constant  which  can  be  evaluated  by  enforcing  (3): 

0  =  D  •  C'  =  D  •  C'  -  d  C'  •  C', 

so  that 


On  combining  (16)  with  (18),  we  obtain 

D  =  D  -  ^  C'  (6  ■  C') 

=  { i  - 1 C'  ®  C’}  [d] 
3  '[*]• 


(16) 


(17) 


(18) 


(19) 


Here  I  is  the  fourth  order  identity  tensor  and  P  is  the  local  fourth  order  “projec¬ 
tion”  tensor  mapping  TJ  into  TLc-  The  geometric  interpretation  of  (16-19)  is  shown 
schematically  in  Fig.  1. 

We  note  in  passing  that  D,  as  decned  by  (19),  is  the  unique  element  of  closest 
to  D  (min  |D  -  D|).  For  those  crystals  having  D  •  C'  =  0,  D  =  D,  analogous  to  the 
TAYLOR  assumption  in  unconstrained  crystals. 
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Now,  by  inserting  (19)  into  the  global  average  compatibility  condition  (14),  and 
recalling  that  D  is  constant,  we  obtain 

<  P  [f>]  >  =  <  P  >  [6]  =  D.  (20) 

A  comparison  of  (20)  with  (15)  indicates  that 

K  =<  P  >~\  (21) 

so  that  local  deformation  is 

D  =  [P  <  P  >-'}  [D]  .  (22) 

We  complete  the  description  of  local  kinematics  by  choosing  the  local  spin  W.  As 
in  TAYLOR  models  which  take  no  account  of  morphological  grain  texture,  we  simply 
equate  local  spin  with  the  global  average  spin: 

W  *  W,  (23) 

trivially  satisfying  <  W  >  =  W.  The  rate  of  change  of  lattice  orientation  is  given  by 
a  spin  tensor  W  -  Wp  which  provides  the  following  evolution  equations  (ASARO  and 
RICE,  1977)  for  R®  and  for  A®,  the  latter  of  which  is  the  skew  part  of  the  slip  system 
dyadic: 


R*  =  (W  -  WP)R®  -  R“(W  -  Wp), 

(24) 

A“  =  (W  -  WP)A®  -  A®(W  -  Wp), 

(25) 

where  the  plastic  spin,  Wp,  within  a  crystal  is  given  by 

it  1 

Wp  =  £f,i{na®s<’-sa®na} 

a=l  2 

=  £if“  A®.  (26) 

ar=l 
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We  note  in  passing  that  the  tensor  C',  being  connected  to  the  crystal  lattice,  also 
evolves  with  this  spin: 


C' 


(W  -  WP)C'  -  C'(W  -  Wp). 


(27) 


3.2  Stress  and  Equilibrium 

The  macroscopic  stress  deviator  S  must  also  be  the  volume  average  of  the  local 
stress  deviators: 

<  S  >  =  S.  (28) 

On  inserting  the  decomposition  for  S  given  by  (8)  into  (28),  we  obtain 

<  S*  >  <  (S  •  C')C'  >  =  S.  (29) 

The  local  stresses  S*  are  determined  in  terms  of  D  by  inserting  (22)  into  (ll)  and 
then  inverting  the  latter  equation.  However,  the  local  deviatoric  stress  component 
Sc  =  |  SC'  is  a  reaction  enforcing  the  inextensibility  constraint  and  can  be  determined 
only  from  equilibrium  considerations.  We  approximate  this  local  stress  component  as 
the  corresponding  projection  (component)  of  the  macroscopic  stress,  S: 

Sc  =  l  S  C'  =  ^S  C\  (30) 

When  (30)  is  inserted  into  (29),  and  the  equation  is  rearranged  (bearing  in  mind  that 
S  is  constant),  there  obtains 

<s->  =  {z-?  <C'®C'>}[§! 

=  <  P  >  (SJ.  (31) 

Equation  (31)  provides  the  connection  between  global  deformation  rate  and  stress 
deviator.  As  has  been  emphasized  by  others  (e.?.,  ASARO,  1983a, b),  the  simulation  of 
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deformation  and  texture  evolution  in  general  should  not  proceed  by  prescribing  a  priori 
all  components  of  D.  Rather,  specific  components  of  D  and  of  S  must  be  given  in  a 
manner  consistent  with  a  well  posed  boundary  value  problem.  For  example,  in  simple 
compression,  one  might  prescribe  £)3S  =  E  ^  0  and  require  Ssi  —  S3  2  —  Sn  —  S22  = 
S12  =  0.  In  any  event,  the  five  prescribed  macroscopic  components,  in  conjunction 
with  the  five  relations  (31),  provide  the  means  for  determining  all  stress  deviators  and 
deformation  rates,  both  locally  and  in  global  average. 

S.S  Tensorial  Representation  of  Texture 

In  our  constrained  hybrid  model,  the  macroscopic  texture  tensor  <  P  >  plays  a 
prominent  role  in  defining  both  the  local  kinematics  and  global  equilibrium.  This 
tensor  can  be  explicitly  evaluated  by  replacing  the  volume  integral  in  its  definition 
by  a  weighted  integration  over  orientations  defined  on  the  unit  sphere,  S.  ON  AT  and 
LECKIE  (1988)  have  recently  given  a  convenient  tensorial  representation  for  a  scalar 
field  on  S.  Here  we  adapt  those  results  for  the  polycrystalline  aggregate. 


Let  an  ensemble  of  Na  grains  of  equal  volume  constituting  a  material  point  have 
c-axis  unit  vectors  (±Ci,  ±C2,  . . . ,  ±c^),  where  we  note  that  the  ±  sign  is  necessary 
since  the  symmetry  of  an  axis  emphasizes  that  no  unique  sign  can  be  assigned  to  c<. 
The  endpoints  of  this  set  of  unit  vectors  pierce  S  at  2 Ng  locations.  We  describe  the 
cloud  of  points  on  S  by  an  axis  orientation  density  function,  i/c(n).  Consider  a  given 
point,  n,  on  S  and  surround  it  with  an  infinitesimal  area  dA.  Let  d Nt  be  the  number 
of  grains  with  c-axis  projecting  into  d A;  the  point  density  dN9/AA  is  normalized  to 
obtain  t'c(n): 


^c(n)  = 


2tt  d N9 
N~g~6A' 


(32) 
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ON  AT  and  LECKIE  (1988)  have  shown  that  functions  uc  of  this  sort  can  be  repre¬ 
sented  in  the  form 

i'c(n)  =  1  +  n)  +  £iju%ju(n)  +  . . .  (33) 

where 

C[- :  cartesian  components  of  a  traceless,  symmetric 
second  order  tensor  (  C,'y  =  e<  •  C'e;) 

Cijki-  cartesian  components  of  a  fourth  order 

completely  symmetric,  completely  traceless  tensor, 
etc., 

and  the  even  Fourier  “basis”  terms  are  modified  spherical  functions  having  cartesian 
components 

(®)  =  ~  2  &ij >  (34) 


Tijki{n)  ~  < 


mn}nkni 

{&ijnkni  +  Suriiitj  +  6iknjrii  +  ^n,nfc  +  6>Jkn,n,  +  ^/n,n*) 
+  5x7  +  Sik&jl  +  tiilfyk) 


(35) 


etc. 


Here  6,,  is  the  Kronecker  delta,  and  n,  are  cartesian  components  of  the  unit  vector 
n  relative  to  the  orthonormal  basis  e,.  In  direct  tensor  notation,  for  example,  F(n)  = 
n  ®  n  -  4 1. 

The  basis  tensors  F,  7,  . . .  generate  their  corresponding  Fourier  coefficients  through 
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Mn): 


1  =  £/,**(»)  X  KM. 

(36) 

C*  =  4?r  x  2  /s^c(n)  x  Fi,(n)dA, 

(37) 

t  1  3  x  5  x  7  x  9  f  .  ,  .  ,  . 

=  4?r  X  2X3X4 

(38) 

etc. 

The  Fourier  coefficients  C/y,  Ci,u,  •  • .  can  be  thought  of  as  “moments”  of  the  axis  orien¬ 
tation  distribution  function,  t»c(n).  ONAT  and  LECKIE  emphasize  that  the  tensor  char¬ 
acter  of  the  Fourier  coefficients  is  such  that,  under  orthogonal  transformations  Q,  they 

_  t  _  t 

transform  via  standard  tensor  transformations  of  appropriate  order,  e.g.,  C  — *•  QC'Q  , 
etc.  For  isotropic  distributions  of  the  c-axis,  v e(n)  =  1,  and  all  other  texture  tensors 
(O',  C,  etc.)  vanish. 

We  can  use  the  Fourier  tensors  to  calculate  those  ensemble-average  tensor  functions 
of  the  texture  which  depend  only  on  the  local  c-axis  orientation  by  evaluating  the 
integral  over  5  of  the  function,  weighted  by  i/c(n);  i.e., 

<  f(c)  >  =  j-  f  i/c(n)  f(n)  dA  (39) 

4tt  Js 

for  any  tensor-valued  function  r(c).  In  the  particular  case  of  the  projection  tensor  P, 

P(C)  =  1-1  (c®c-il)®(c®c-jl) 

=  J-?P(c)«F(c).  (40) 

We  can  define  the  average  value,  <  P  >,  by  its  action  on  an  arbitrary  G  6  Ij.  On 
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using  (33-40),  the  result  can  be  expressed  as 


<  P  >  [G]  = 


I  3x4 
r 5x7x9 


(C'G)I— sMc'G  +  GC'} 


(41) 


3x4 

5x7x9 


C  [G 


precisely,  independent  of  any  higher  Fourier  moment  of  texture.  For  isotropic  distribu¬ 
tion  of  the  c-axis,  C*  =  C  =  0.  In  such  cases,  for  example,  from  (15),  (21),  and  (41), 
<  P  >  [D]  =  |  D  so  that  D  =  |  D. 


The  tensor  mapping  on  the  unit  sphere  presumes  a  smoothed,  infinite  dimensional 
basis  for  the  axis  orientation  distribution,  whereas  the  motivation  (and,  indeed,  nu¬ 
merical  simulations  of  polycrystals)  are  based  on  a  finite  number  of  grains.  Rational, 
effective  means  for  obtaining  the  smoothed  texture  tensors  C',  C,  etc.,  from  finite  pop¬ 
ulations  have  been  developed  based  on  invariant  tesselation  and  quadrature  on  the 
manifold  S,  PARKS  (1989). 


8.4  Generalization  to  Two  or  More  Constraints 

The  method  outlined  in  sections  3.1  and  3.2  is  applicable  to  polycrystalline  models 
of  crystals  having  a  single  kinematic  deficiency  in  deviatoric  response.  We  briefly  gen¬ 
eralize  the  formulation  to  accommodate  crystals  having  two  linearly  independent  con¬ 
straints  on  deviatoric  deformation;  further  generalization  to  three  or  more  constraints 
is  transparent.  The  extended  model  forms  the  basis  for  simulations  of  deformation  and 
texture  evolution  in  polycrystalline  olivine  reported  in  section  4.3,  below. 

Suppose  that  deformation  in  a  crystal  satisfies  the  condition  D  •  B'  =  0  in  addition 
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to  (3),  where  the  deviatoric  tensor  B'  is  orthogonal  to  C':  B'  •  C'  =  0.  We  again 
introduce  the  intermediate  tensor  D  using  (15),  but  generalize  (16)  to 

D  =  D-dC'  —  6B'.  (42) 


Enforcement  of  the  two  constraints  provides  b  =  (D  •  B')  / (B;  •  B')  as  well  as  the 
result  (18)  for  d.  The  generalization  of  (19)  is 

D  =  £[D],  (43) 

where 

Satisfaction  of  (14)  now  requires  <  Jl  >-1=  K. 


The  invariant  decomposition  of  stress  deviator  analogous  to  (8)  is 


S  =  S"  +  ScC’  +  SbB' 


(S-C')  (S-B-) 

+  (C'-C')  (B'  ■  B')  ’ 


(45) 


where  the  kinematically-determinable  S**  lies  in  a  three-dimensional  subspace  of  T-J, 
and  Sc  and  Sb  are  the  deviatoric  stress  reactions  which  must  be  determined  from 
equilibrium.  Finally,  the  approximations  S  •  C'  =  §  •  C;  and  S  •  B'  =  S  •  B'  furnish, 
through  (28),  the  generalized  global  equilibrium  relation 


<  Z  >  (S)  =<  s**  >  . 


(46) 
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4.  Applications 


4.1  Numerical  calculations 

Because  of  the  incompressibility,  the  strain  rate  is  vectorized  to  five  independent 
components  referred  to  a  fixed  cartesian  basis: 

D  — »  D  —  |(D2j  —  Du)/ \/2,  y^3/2jD3j,  jDjj/ v^2,  Du/ \/2,  D«/ \/2  j  . 

All  other  deviatoric  tensors  are  also  vectorized  in  the  above  convention. 

The  solution  procedure  is  as  follows.  Once  an  estimate  of  ©  is  known,  the  local 
stretching  D  is  calculated  from  (22),  since  both  local  orientation  and  global  texture  are 
assumed  to  be  known  instantaneously.  The  non-linear  relation  (11)  between  D  and  the 
reduced  stress  S*  is  solved  by  the  Newton-Rap hson  method  to  obtain  S*.  The  global 
average  <  S*  >  is  then  obtained,  and  (31)  is  solved  to  obtain  S.  If  global  stress  bound¬ 
ary  conditions  are  satisfied  to  sufficient  tolerance,  the  solution  is  accepted;  otherwise, 
Newton-Raphson  corrections  to  the  work-conjugate  components  of  D  are  made,  and 
the  procedure  is  repeated.  Once  global  conditions  are  met,  the  remaining  portion  of  S  is 
obtained  from  assumption  (30)  when  it  is  required.  Local  updating  of  lattice  orientation 
proceeds  by  constructing  an  orthogonal  tensor  Q1,  =  exp  {At  (W  -  Wp)},  where  At  is 
a  small  time  increment,  and  the  spin  tensors  are  obtained  from  (23)  and  (26).  Lattice 
tensors  ( e.g Ra,  A®,  C')  are  then  updated  according  to  C'(t  -f  At)  —  QLC'(t)QLT , 
etc.  In  the  updated  configuration,  the  procedure  is  repeated  for  the  next  increment. 
Typically,  global  deformation  increments  of  one  percent  are  used.  Similar  procedures 
are  used  for  the  case  of  two  constraints. 
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Initial  isotropic  distributions  of  lattice  orientation  are  constructed  using  random 
numbers.  A  check  on  the  isotropy  of  the  initial  texture  was  obtained  from  (41);  typically 
the  numerical  value  of  <  P  >  deviated  from  f  J  by  less  than  two  percent. 

4-2  Application  to  Hexagonals 

The  slip  systems  reported  to  be  activated  in  different  hexagonal  crystals  (GROVES 
and  KELLY,  1963)  axe  of  two  types  —  slip  systems  with  a  as  Burger’s  vector:  basal 
slip  {0001}  <  2110  >,  prismatic  slip  {01 10}  <  2110  >,  and  pyramidal  slip  {0111} 
<  2110  >;  and  slip  systems  with  (a  +  c)  as  Burger’s  vector.  Tom£  and  KOCKS  (1985) 
analyzed  single  crystal  plasticity  of  hexagonals  yielded  by  different  types  of  slip  and 
twinning  systems. 

Since  the  Burger’s  vector  a  of  basal  and  prismatic  systems  lies  in  the  basal  plane, 
these  two  sets  of  systems  together  comprise  only  four  linearly  independent  slip  systems. 
For  the  same  reason,  the  addition  of  the  pyramidal  slip  {0111}  <  2110  >  does  not  sup¬ 
ply  the  missing  degree  of  freedom  (».«.,  no  plastic  deformation  along  the  c-axis).  The 
additional  mechanisms  which  can  supply  the  missing  degree  of  freedom  are  pyramidal 
slip  with  Burger’s  vector  <  a  -I-  c  >  or  certain  twinning  systems,  TOME  and  KOCKS 
(1985). 

HUTCHINSON  (1977)  applied  self-consistent  and  upper  bound  models  to  hexagonal 
polycrystals  with  isotropic  texture  and  with  an  ideal  lattice  ratio  c/a  =  2y(2/3)  to 
predict  the  uniaxial  reference  stress  oQ  in  the  following  uniaxial  power  law  creep  equa¬ 
tion: 
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D"  =  io(oe'/e0)n, 


(47) 


or,  equivalently, 

d0  =  oe'(lr'lio)-lln,  (48) 

where  Deq  =  D  •  D  and  =  ^/|S  •  S  are  the  macroscopic  equivalent  uniaxial  ten¬ 
sile  strain  rate  and  equivalent  tensile  stress,  respectively,  70  is  the  reference  shear  rate 
of  the  single  crystal,  and  n  =  1/m,  where  m  is  the  rate  sensitivity  exponent.  HUTCHIN¬ 
SON  considered  basal  and  prismatic,  as  well  as  pyramidal  slip  {1122}  <  1123  >.  r  ' 
compared  self-consistent  and  upper  bound  (TAYLOR)  results  for  different  rate  sensit 
ities  and  differing  ratios  of  critical  resolved  shear  stresses  on  different  systems.  When 
the  critical  resolved  shear  stress  on  the  pyramidal  system  is  infinite  (t'.e.,  pyramidal 
slip  is  omitted),  the  upper  bound  calculation  fails  because  of  the  inextensibility  con¬ 
straint  in  the  c  direction.  However,  in  this  case  he  was  able  to  obtain  results  from  the 
self-consistent  model. 

Here  we  compare  results  from  the  CH  model  with  HUTCHINSON’S  in  the  case  of 
four  independent  systems  only  {basal  and  prismatic  systems).  Direct  comparisons 
we  made  for  exponents  n  =  1  and  n  =  3  and  for  the  critical  resolved  shear  stress 
ratio  Tpjrh  =  10,  where  r*  and  rp  are  critical  resolved  stresses  for  basal  and  prismatic 
systems,  respectively.  The  polycrystalline  aggregate  consists  of  248  randomly  oriented 
single  crystals,  resulting  in  an  isotropic  texture. 

Results  of  the  tensile  reference  stress  d0,  normalized  by  rp,  are  given  in  Table  1, 
as  are  HUTCHINSON’S  results.  In  the  linear  (n  =  1)  case,  the  model  prediction  falls 
between  the  self-consistent  and  Hashin-Shtrikman  lower  bound,  approximately  78%  of 
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the  self-consistent  value.  For  n  =  3,  the  self-consistent  estimate  more  than  doubles, 
compared  to  the  n  =  1  case,  while  the  CH  estimate  increases  modestly  by  only  ten 
percent  above  its  n  =  1  value.  HUTCHINSON  (1977)  noted  that  there  are  reasons  to 
suspect  that  the  self-consistent  model  overestimates  o0  when  pyramidal  slip  is  absent, 
so  the  current  results  may  be  more  accurate. 

Figure  2  shows  the  variation  of  the  normalized  reference  stress  versus  the  slip  system 
hardness  ratio  r*/rp  for  n  =  1  and  n  =  3  as  predicted  by  the  present  model.  Results 
for  n  >  3  differ  only  slightly  from  those  of  n  =  3.  We  note  that  the  value  of  n  =  3 
and  “hard”  prismatic  systems  (small  r4/rp)  corresponds  to  ice  crystals  (HUTCHINSON, 
1977;  WEERTMAN,  1973;  ASHBY  and  DUVAL,  1985). 

4-8  Application  to  Olivine 

Olivine  is  a  crystalline  mineral  which  is  believed  to  be  the  major  constituent  of 
the  upper  mantle.  Strong  seismic  anisotropy  has  been  observed  in  the  upper  mantle 
beneath  both  oceanic  and  continental  crust.  At  high  pressures,  where  microcracks 
are  closed,  it  is  possible  that  seismic  anisotropy  could  be  due  to  preferred  crystalline 
orientation,  or  texture.  Recently  TAKESHITA,  et  al.,  (1989)  have  begun  to  model  the 
development  of  polycrystalline  deformation  texture  in  olivine. 

Olivine  has  an  orthorhombic  crystal  lattice.  Several  experiments  have  been  carried 
out  to  determine  operative  slip  systems  at  different  temperatures  (CARTER  and  Ave 
Lallemant,  1970;  Durham  and  Goetz,  1977;  Raleigh,  1968).  Takeshita,  et  al., 
(1989)  summarize  the  slip  systems  into  two  classes;  class  A  for  low  temperature,  and 
class  B  for  higher  temperature,  and  estimate  the  resistances  for  these  systems.  They 
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note  that  the  class  B  systems  and  resistances  are  based  on  single  crystal  experiments, 
while  the  class  A  systems  should  be  viewed  as  more  tentative.  These  systems  are 
summarized  in  Table  2.  Both  classes  have  only  three  independent  slip  systems,  allowing 
no  inelastic  stretching  along  (unit)  crystallographic  directions  a,  b,  and  c  (».e..,  normal 
strain  components  referred  to  crystallographic  axes  cannot  be  accommodated  by  slip). 

In  order  to  use  the  TAYLOR  formalism,  TaKESHITA,  et  al.,  proposed  that  the  normal 
strains  (strains  in  crystallographic  directions)  are  accommodated  by  dislocation  climb 
while  shear  strains  (as  referred  to  crystallographic  axes)  axe  accommodated  by  slip.1 
The  neglected  normal  strain  components  were  presumed  to  make  no  contribution  to 
textural  evolution.  They  also  used  the  self-consistent  model  of  MoLINARI,  et  al., 
(1987)  to  predict  texture  of  olivine  deformed  by  simple  shear  and  compared  results  to 
the  TAYLOR  model  predictions.  For  neither  model  was  the  predicted  polycrystalline 
stress  strain  response  reported. 

Here  we  use  the  nonhardening  CH  model  to  predict  the  behavior  of  olivine  and 
compare  our  results  to  those  of  TAKESHITA,  et  al.  Since  there  are  only  three  indepen¬ 
dent  slip  systems  for  either  of  the  two  classes  of  slip  systems,  we  use  the  two-constraint 
version  of  the  model  outlined  in  section  3.4.  The  normal  stretching  components  are 
zero  along  the  crystallographic  axes:  a  ■  Da  =  b  •  Db  =  c  •  Dc  =  0.  Because  of  as¬ 
sumed  incompressibility,  the  three  normal  strain  constraints  are  reduced  to  only  two 
independent  deviatoric  constraints,  which  c?n  be  expressed  as 

D  •  B'  =  0,  (49) 

'Procedures  locally  requiring  certain  deformation  components  to  equal  the  corresponding  macroscopic 
average  components,  ■while  simply  igm  -ing  (in  a  mechanics  sense)  the  other  components,  have  collectively 
come  to  be  termed  “relaxed  constrai.  models. 
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D  •  C'  =  0, 


(50) 


where  B'  =  (b  ®  b  —  a  ®  a)  and  C'  is  given  by  (3). 

We  considered  the  olivine  behavior  and  development  of  anisotropy  under  conditions 
of  both  simple  shearing  and  plane  strain  compression  at  high  temperature  (class  B). 
The  initial  texture,  consisting  of  248  randomly  oriented  single  crystals,  is  isotropic  and 
is  represented  by  pole  figures  in  t  ig.  3.  We  chose  a  value  of  n  =  5  as  did  TAKESHITA, 
et  al.,  (1989).  The  predicted  stress  strain  curves  are  given  in  Fig.  4  for  both  plane 
strain  compression  and  simple  shear.  The  macroscopic  equivalent  tensile  stress,  de?, 
is  normalized  by  r0,  the  initial  resistance  of  the  easiest  slip  system.  The  macroscopic 
equivalent  strain  is  defined  as  c*4  =  /  D eq  df . 

The  texture  developed  under  simple  shear  is  shown  in  Fig.  5,  where  the  sense  of 
shear  is  indicated  by  the  arrows.  This  texture  is  in  agreement  with  those  obtained 
by  TAKESHITA,  et  al.,  and  with  a  texture  which  has  been  observed  in  naturally  de¬ 
formed  olivine  (MERCIER,  1985).  Texture  developed  during  plane  strain  compression, 
Fig.  6,  as  predicted  by  the  present  model,  is  also  in  agreement  with  the  predictions  of 
TAKESHITA,  et  al.  Unfortunately  for  purposes  of  comparison,  at  present  there  seems  to 
be  no  available  experiments  on  polycrystalline  olivine  documenting  texture  evolution 
from  an  initially  isotropic  texture  under  simple  shear  or  plane  strain  compression. 

4-4  Application  to  HDPE 

High  density  polyethylene  (HDPE)  is  a  semi-crystalline  polymer  with  crystallinity 
typically  between  70%  and  80%.  In  order  to  apply  the  present  model  in  a  highly 
idealized  manner  to  HDPE,  we  neglect  the  contribution  of  the  amorphous  phase  by 
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assuming  an  “ideally  crystalline”  polyethylene  (100%  crystallinity).  We  further  neglect 
the  important  morphological  texture  of  the  material,  which,  in  the  macroscopically 
isotropic  case,  usually  consists  of  spherulitic  regions  containing  crystalline  lamellae  and 
amorphous  interlamellar  domains  of  highly  correlated  local  orientation,  but  of  random 
orientaion  when  averaged  over  the  spherulite. 

Plastic  deformation  of  single  crystal  (lamella)  HDPE  was  reviewed  by  BOWDEN  and 
YOUNG  (1974)  and  by  HAUDIN  (1982).  The  crystalline  domains  (lamellae)  form  an 
orthorhombic  lattice  with  the  c  axis  in  the  polymer  chain  direction.  The  crystals  can 
deform  by  slip,  twinning  or  martensitic  transformation.  In  this  application  we  consider 
only  plasticity  by  slip.  The  possible  slip  systems  are  chain  slip  (axis  c  is  Burger’s  vector) 
and  transverse  slip  (Burger’s  vector  perpendicular  to  c  axis).  Table  3  summarizes 
all  possible  slip  systems,  along  with  corresponding  estimated  normalized  deformation 
r»*istances.  Together,  these  provide  only  four  independent  slip  systems,  allowing  no 
plastic  extension  in  the  c  direction  (chain  direction).  The  possible  twinning  modes, 
which  are  not  included  in  our  calculations,  do  not  provide  an  additional  independent 
system. 

We  used  the  CH  model  with  one  constraint  to  predict  the  behavior  and  texture 
evolution  of  an  ideally  100%  crystalline  HDPE  under  both  tension  and  simple  shear. 
The  initial  texture,  again  consisting  of  248  randomly  oriented  single  crystal  orientaions, 
is  isotropic,  as  indicated  by  the  pole  figure  shown  in  Fig.  3.  The  value  of  n  is  taken  as 
17,  which  corresponds  approximately  to  the  strain  rate  sensitivity  exponent  m  =  0.06 
obtained  experimentally  at  room  temperature  by  G’SELL  and  JONAS  (1979)  in  tests 
of  spherulitic  HDPE.  In  view  of  the  small  lamellar  thicknesses  (~  20 nm),  it  was  felt 
that  strain  hardening  by  dislocation  interaction  within  the  lamellae  was  likely  to  be 
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insignificant,  so  no  hardening  of  slip  system  critical  resolved  shear  stress  was  included. 


Equivalent  stress  strain  curves  for  simple  shear  and  for  tension  are  given  in  Fig.  7. 
Stress  is  normalized  by  the  critical  resolved  shear  stress  for  chain  slip.  Evidently, 
textural  hardening  of  HDPE  is  very  important  in  tension  but  not  in  simple  shear.  This 
result  is  in  excellent  agreement  with  the  experimental  results  of  G  ’SELL,  et  al.,  (1982). 

The  texture  developed  during  tension  is  a  fiber  texture,  Fig.  8,  with  chain  direction 
(c-axis)  parallel  to  the  tension  direction.  Clearly,  this  distribution  of  the  inextensible 
chain  directions  is  responsible  for  the  dramatic  upturn  in  the  tensile  stress  strain  curve 
of  Fig.  7. 

The  texture  developed  by  simple  shear  is  shown  in  Fig.  9,  where  it  is  noted  that 
a  reversed  sense  of  shear  (as  compared  to  that  used  for  olivine  in  Fig.  5)  has  been 
adopted.  In  Fig.  9a,  at  a  shear  of  7  =  2.16,  a  sharp  c-axis  texture  component  forms 
along  the  material  fiber  of  maximum  stretch,  and  the  two  other  crystallographic  axes  are 
almost  randomly  oriented  around  c  (fiber  texture).  Figure  9b,  at  a  shear  of  7  =  2.6, 
shows  further  sharpening  of  the  c-axis  texture  in  a  direction  closer  to  the  shearing 
direction.  The  fact  that  the  c-axis  rotates  toward  the  direction  of  maximum  stretch, 
which  itself  is  tending  toward  the  shearing  direction  with  ongoing  shear,  allows  most 
of  the  deformation  to  be  accomodated  by  chain  slip  (easy  slip).  This  explains  why 
the  textural  hardening  is  very  weak  in  simple  shear.  On  the  contrary,  when  the  HDPE 
becomes  textured  under  tension,  it  exhausts  available  deformation  mechanisms  and  the 
requisite  driving  stress  tends  to  infinity. 

Results  from  more  sophisticated  composite  inclusion  models  (AHZI,  PARKS  and 
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argon,  1989;  PARKS  and  AHZI  1989)  i„  k 

lamellar  amorphous  phase  in  HOPE  eh  '  d  "r  ”  “C°“med  *”  CffeClS  °f  the  inter- 
were  obtained  1  •  '  e<1,,iVate,,t  ‘«.dree 

::  r  c~y  ^ *-  -  *•  «-  — —  ** 

primary  ro!e  of  the  amorphone  materia,  in  HOPE  „  t0  provide 
straining  without  eubetantia!  effect  on  crystallographic  texture. 


5.  Discussion 


We  have  clearly  identified  the  special  constitutive  features  of  rigid-viscoplastic  single 
crystals  which  do  not  possess  five  independent  slip  systems.  In  addition  to  the  rela¬ 
tively  well-known  kinematic  constraint,  there  is  a  corresponding  indeterminacy  in  the 
deviatoric  stress.  Based  on  these  observations,  we  have  constructed  a  simple  TAYLOR- 
type  polycrystalline  model  which  accounts  for  the  local  kinematic  constraints,  as  well 
as  providing  a  means  for  dealing  with  the  “kinematically  indeterminant”  part  of  the 
local  stress  deviator.  As  examples  of  the  broad  applicability  of  the  model,  we  applied 
it  to  various  problems  of  polycrystalline  deformation  and  texturing  in  kinematically 
deficient  crystals  of  hexagonal  and  orthorhombic  symmetry.  These  examples  were  in¬ 
tended  simply  to  demonstrate  the  new  model,  rather  than  as  comprehensive,  or  even 
complete,  studies  of  particular  problem  areas.  Further  applications  of  the  model  are 
underway,  and  will  be  reported. 

A  primary  reason  for  the  success  of  the  model  is  that,  like  the  self-consistent  formu¬ 
lations,  it  imposes  far  less  stringent  kinematic  conditions  on  individual  inclusions  than 
uniform  deformation.  As  in  the  self-consistent  models,  global  compatibility  (in  the 
form  of  macroscopic  deformation  equalling  the  volume  average  of  local  deformation)  is 
preserved,  but  local  compatibility  is  not  addressed.  While  such  “disregard”  for  local 
compatibility  may  seem  to  be  a  limitation  of  these  approximate  models,  a  rigorous 
study  would  require  higher  order  statistical  information  regarding  the  microstructure, 
2is  well  as  recognizing  the  inhomogeneity  of  plastic  deformation  within  grains  (ASARO, 
1983a, b). 

In  preserving  the  volume-average  character  of  deformation,  a  macroscopic  projec- 


tion  tensor,  <  P  >,  played  a  prominent  role.  We  presumed  that  this  diagonally  sym¬ 
metric  fourth  order  tensor  was  invertible;  ».e.,  that  it  mapped  Tj  to  itself.  In  certain 
extreme  (Dirac-like)  textures,  such  as  the  fiber  textures  displayed  in  Figs.  8  and  9, 
<  ?  >  itself  becomes  degenerate,  spanning  only  TjXC.  The  simple  shear  calculations 
could  be  continued  without  substantial  difficulty,  but  tensile  deformation  was  extin¬ 
guished.  In  any  event,  the  numerical  condition  of  <  P  >  can  become  an  issue  at  sharp 
texture. 

The  tensor  <  P  >  also  played  an  important  role  in  enforcing  the  volume  average 
condition  on  the  stress  deviator.  This  occurred  because  the  local  constraint  stress  was 
approximated  by  the  corresponding  component  of  the  macroscopic  stress.  We  are  not 
aware  of  rigorous  justifications  for  making  this  assumption.  Rather,  we  sought  simply 
to  provide  for  some  measure  of  equilibrium  between  the  constrained  inclusion  and  its 
macroscopic  surroundings.  Perhaps  a  related  (though  less  sophisticated)  assumption 
might  attempt  to  account  for  pressure  sensitivity  of  rigid-viscoplastic  polycrystalline 
slip  by  assuming  local  pressure  within  an  inclusion  to  equal  the  global  average  pressure, 
which  in  turn  could  be  deduced  from  traction  boundary  conditions. 

We  generalized  the  formulation  to  account  for  two  (or  more)  deviatoric  constraints 
within  grains  (three  or  fewer  independent  systems),  and  applied  it  to  deformation 
texturing  in  the  mineral  olivine.  Formally,  the  methodology  could  be  extended  to  a 
single  slip  system  within  each  grain,  but  we  reserve  judgment  on  its  utility  in  that 
limit.  For  example,  in  the  hexagonal  reference  stress  study,  we  found  indications  of 
a  non-zero  limit  for  tf0/rp  48  r*/rj>  approached  zero,  Fig.  2.  On  the  other  hand,  one 
could  argue  that  in  this  limit,  prismatic  slip  would  be  “infinitely”  more  difficult  than 
basal  slip,  so  that  only  the  two  independent  basal  systems  should  be  retained  (three 
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constraints).  In  that  event,  90  would  scale  with  the  retained  rt. 
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(rt/rp  =  1/10) 

oo hv  (n  =  !) 

— 

/rp  (n  =  3) 

Taylor  Model 

oo 

oo 

H.-S.  Upper  Bound 

00 

oo 

Self-Consistent 

1.75 

3.9 

H.-S.  Lower  Bound 

0.74 

Sachs  Model 

0.45 

CH  Model 

1.37 

1.49 

Table  1:  Macroscopic  tensile  stress  reference  values,  normalized  by  rp, 
in  an  isotropic  polycrystal  of  hexagonal  single  crystals  without 
pyramidal  slip.  Comparison  of  CH  predictions  to  results  of 
HUTCHINSON  (1977)  for  the  case  rt/rp  =  1/10. 


slip  system 

normalized  resistance 

class  A 

{110}[001] 

1. 

{011}[100] 

1. 

(010)  [100] 

1. 

class  B 

(001)  [100] 

1.3 

(010)  [001] 

2.7 

Table  2:  Slip  systems  of  olivine  and  their  corresponding 
normalized  initial  critical  resolved  shear  stresses. 

A:  Low  temperature;  B:  High  temperature 
(after  TAKESHITA,  tt  a/.,  1989). 
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slip  system 

normalized  resistance 

(100)[001) 

1. 

Chain  slip 

(010)[00l] 

1. 

{110}[001] 

1. 

(100)  [010] 

1.1 

Transverse  slip 

(010)[100] 

1.4 

{110}  <  110  > 

1.8 

Table  3:  Slip  systems  of  polyethylene  and  estimates  of  their  corresponding 
normalized  initial  critical  resolved  shear  stresses. 


Figure  2:  Macroscopic  reference  stress  for  isotropic 
HCP  polycrystals  having  only  basal  and  prismatic  slip. 


(002) 


Figure  4:  Macroscopic  equivalent  stress-strain 
curves  for  olivine. 


(001) 


(010) 


Figure  5:  Predicted  (001)  and  (010)  pole  figures 
for  simple  shear  of  olivine  (y=1.3). 
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